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Abstract 

The MINOS experiment has observed a rise in the underground muon 
charge ratio = /i'^/yU". This ratio can be related to the atmospheric 
production ratios of tt^/tt^ and K~^/K~. Our analysis indicates that the 
relevant variable for studying the charge ratio is E^^'^-l'"''^'^ cosd, rather than 
j^surface^ We Compare a simple energy dependent parameterization of the 
rise in the charge ratio with more detailed previously published Monte Carlo 
simulations and an analytical calculation. We also estimate the size of two 
previously neglected effects in this context: the charge sign dependency of 
the dE/dx in rock, and the energy dependence of heavy primaries on the 
derived K~^/K~ ratio. 
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1. Importance of Charge Ratio Measurements 

Atmospheric muons come dominantly from the decay of vrs and Ks pro- 
duced in hadronic showers when cosmic rays interact in the earth's atmo- 
sphere. These muons have been studied with energies ranging from hundreds 
of MeV to well over a TeV. A quantitative understanding of cosmic ray muons 
has value for a number of diverse topics, from atmospheric neutrinos to the 
chemical composition of the highest energy cosmic rays. The charge ratio of 
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cosmic ray muons has been previously measured over three orders of magni- 
tude in energy. Recently the MINOS experiment [l|, ^ presented data that 
for the first time showed a rise in the measured charge ratio 



from previous measurements at high values of or more specifically, high 
values of ii;^"''/'^ce ^ 

In this paper, we discuss some of the issues involved in the measurement 
and interpretation of the muon charge ratio. In particular, we develop a 
simplified model where the rise in the charge ratio can be understood from the 
properties of vr and K mesons, and the observation of the rise can be used to 
determine the / vr" and / K~ ratios. We also address several other issues 
related to the measurement of the charge ratio, including the role of muon 
energy loss, a detector's Maximum Detectable Momentum (MDM), and the 
effect of possible differences in the spectral index of cosmic ray Hydrogen and 
Helium on the interpretation. 

Since the primary cosmic rays are mostly positively charged protons, more 
secondary 7r+ are expected than vr". The quark content of the protons and 
of the atmosphere has been used to estimate the tt^/tt" ratio to be near 
1.27 j^. The charge ratio for kaons is even higher due to the phenomenon of 
associated production. Strange particle production starts with the creation 
of an s quark and an s quark. An s quark which ends up in a nucleus 
is associated with an s quark in a K^{su). The s quark will not be in a 
baryon. There is also K^K~ pair production. Phase space favors hadronic 
production of K^K over pairs at all energies, so large /K~ ratios 

are expected. 

A standard parametrization of the atmospheric muon energy spectrum is 
given by Gaisser[3]: 



The two terms inside the curly bracket give the contributions from charged 
pions and kaons. For the hadron mass {rrii) and lifetime (r^), and the atmo- 
spheric scale height (h ~ 6.4 km), the values of = mich/ri are the energies 
where the probability of meson interaction in the atmosphere and decay are 






The zenith angle is denoted by 
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6. In those experiments with a flat overburden, the largest muon intensity 
is at cos{6) = 1. Reference [3| gives r] = 0.054, which sets the relative con- 
tribution to the muon flux from vr and K decay. The parameter rj depends 
upon the n/K ratio, branching ratios, and kinematic factors which arise due 
to differences between the vr and K masses. 

The ratio between the two terms in Equation [2] quantifies the relative 
contribution to the muon flux from muons due to vr and K decay. It depends 
on £'^"''/"ce cos^ and rises from r] = 5.4% at low energies to rj x ex / = 40% 
at high E'l^^^'^'^^ co^O. This substantial rise which occurs mostly between 115 
GeV and 850 GeV is caused by the different masses and lifetimes of pions 
and kaons, not by any increased amount of kaon production from primary 
cosmic ray interactions in the atmosphere. 

MINOS has provided the first high statistics measurements of the muon 
charge ratio at E^^^^^"-'^^ cos6 > 115 GeV. We will show that this gives the 
needed sensitivity for extracting information about cosmic ray produced pi- 
ons and kaons separately. We also note that a similar calculation for neutrinos 
shows that kaons are the dominant parent for TeV neutrinos, which are the 
largest backgrounds for astrophysical source searches at neutrino telescopes. 

The organization of this paper is as follows: In the following section, we 
review several measurements of the atmospheric muon charge ratio. In Sec- 
tion [31 we discuss some of the issues involved in measuring the charge ratio 
in a magnetic detector underground. In Section HJ we address particular 
issues related to the muon energy loss and how the muon energy at an un- 
derground detector is related to the energy in the atmosphere. In Section |5l 
we use Equation [21 which was based on ideas developed by Zatsepin[5], to 
parameterize the charge ratio in a simplified model. Zatsepin's work is the 
earliest reference that we have found which defines the roles of the criti- 
cal energies = 115 GeV and ex = 850 GeV. There are two important 
implications from this model. First, the charge ratio depends not simply 
on the muon energy i<;^"''/«"=e^ ]^yj^^ ^y^q combination of the energy and 
zenith angle E'^"^-^"^'^ cos^. Second, together with previous measurements at 
low E^^^-l'"''^'^ cos6, observation of the rise in the charge ratio, which mainly 
occurs between e^^ and ex, can be used to fit the meson production charge 
ratios K'^/K~ and /ti~ . In subsequent sections we discuss how this model 
compares to several full Monte Carlo calculations of the muon charge ratio, 
and how a possible difference in the spectral index of the primary cosmic 
Hydrogen and Helium flux might affect the interpretations. 
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2. Compilation of Measurements of the Atmospheric Muon Charge 
Ratio 



Numerous measurements of the atmospheric muon charge ratio data have 
been made at the earth's surface, from MeV to muhi-TeV energies. The 
Hebbeker and Timmermans 2002 compilation article 0] provides many ref- 
erences to these data. Figure 1 displays muon ratio data as a function of the 
surface energy; data from six experiments in the energy range 0.10 TeV to 
10 TeV are shown. 

Baxendale et al. jsj published extensive muon flux data in the momentum 
range 7-500 GeV/c using the Durham spectrometer. This detector consisted 
of magnetized iron blocks, flash tubes, and scintillation counters. The mag- 
netic fleld in the three magnetized blocks was reversed regularly to reduce 
bias effects. Their highest energy charge ratio data point is at 358 GeV/c. 
The publication did not give the ratio as a function of zenith angle, and thus 
only the energy dependency of the charge ratio is used here. 

The CosmoALEPH experiment [7:] at the LEP facility at CERN published 
8 data points in the 80-1600 GeV surface energy range. The cosmic ray 
portion of the experiment used only the hadron calorimeter and the time 
projection chambers. The detector was located at a depth of 320 m.w.e. 
Apparently the magnetic field was not reversed during data collection. The 
zenith range of 0-10 deg was studied, but the publication did not give the 
ratio as a function of zenith angle. Thus only the energy dependency of the 
charge ratio is used here. 

The L3-I-C experiment!^ has published extensive muon charge ratio data 
at many zenith angles at momentum up to 380 GeV/c. The high accuracy 
and extensive angular range of this data are ideal for comparison with our pa- 
rameterization of the muon charge ratio. We have used this data extensively 
in Section 4 of this paper. 



Matsuno et al. [10| published data in the momentum range from 70 
GeV/c to 20 TeV/c using the MUTRON magnetic cosmic ray spectrometer. 
It consists of a solid iron magnetic spectrometer, a calorimeter of 16 pro- 
portional chambers, and 48 spark chambers. Data were collected with both 
forward and reversed magnetic fields. The data were collected only at large 
zenith angles 86-89 degrees. While the data sample is large, there are only a 
few hundred events above 2 TeV/c. 



MINOS has published data [l|, from both their Far Detector at Soudan, 
Minnesota and their Near Detector at Fermilab. Data from both detectors 
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were collected with both forward and reversed magnetic fields. Data are 
available as a function of the zenith angle. 



Rastin published data in the muon surface energy range from 6 - 
1288 GeV using the Nottingham solid iron magnetic spectrometer. Data 
were collected with both forward and reversed magnetic fields. 

Experiments with data points with relatively large error bars have not 
been included in this analysis since they are not useful for the precise studies 
in this paper. Experiments with charge ratio data values at energies above 
100 GeV, with all data points having quoted errors greater than +-0.1, have 
not been included in our analysis; thus Hayman and Wolfendale Ap- 
plet on et al. 13|, Nandi et a l. [l4l| . and Kremer et al. 15| are not included. 
Specifically, Burnett et al. [l6| was not included in our study due to the 
absence of any published estimates of systematic error on their ratio values. 

As seen in Figure [H the ratio has an energy dependency; it is near 1.25 
near 100 GeV and rises to a value near 1.4 at several TeV. The rise in the 
ratio has been attributed to an increasing contribution from kaon decay [3|. 
The data in Figured] could be parameterized with a linear function of E^^'^^"''^^ 
or a log(ii^^"^'^'^'^^) function, but there is no physics reason to expect a linear 
dependency with E, or even a log(E) dependency. We have chosen to use a 
different parameterization that encompasses more of the kinematics and the 
physics of the muon production processes, as will be described in Section [51 



3. Issues in Measuring the Charge Ratio 

In this section, we review the method used for determining the muon 
charge, and consider some of the systematic errors which affect such a mea- 
surement. We describe the concept of Maximum Detectable Momentum, and 
consider two distinct types of systematic errors: bias and randomization. 

Most measurements of the atmospheric muon charge ratio involve the 
use of a magnet to deflect the trajectory of a charged particle. The resulting 
curvature is used to measure both the charge sign and momentum of the 
muon. Depending on the size, granularity, and strength of the magnetic 
fleld, every detector has a Maximum Detectable Momentum (MDM). We 
deflne the MDM as that momentum for which a nearly straight real track 
will have a measured curvature (determined from a flt to points along the 
track) which is one standard deviation from zero. For a uniform magnetic 
fleld, the MDM is simply the reciprocal of the error (s.d.) of the curvature 
measurement, when the curvature is expressed in (GeV/c)~^. 
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Figure 1: Compilation of muon charge ratio from experiments with at least one data point 
with an uncertainty less than 0.10. 



The MDM can be calculated for each event in a detector. As a simplified 
but instructive example, we consider a simplified toroid detector which is 
a long right circular cylinder of radius R, with uniform density, with a B 
field that is azimuthal and constant. The long z-axis is perpendicular to the 
zenith direction. We do not consider any tracks which enter the front or 
leave the rear of the cylinder, i.e. all muons enter and exit the side. All 
such fi tracks will be roughly S shaped, with curvature one way for the first 
half of the track, and the other way for the second half. The length of the 
corresponding half-track is 

L^VR^-bysm{ijj) (3) 

where b is the impact parameter or distance of closest approach to the center 
of the toroid, and ijj is the angle of the track with the z axis (not the zenith 
angle). For this discussion, i/j is near 90 degrees. To calculate the precision 
of the momentum measurement, we need the component of the magnetic 
field B which is perpendicular to the track direction. Approximating at the 
midpoint of the half-track, assuming that ip is not small, we find 

BpERP = 5 X ^[1 - sin2(V;) X cos2(arctan( - &^)/26))] (4) 
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Equation H] is used in Reference [10| , and in Reference [17|] to calculate the 
MDM for MINOS, which varies dramatically between 470 GeV to as low as 
15 GeV for large impact parameters and large ip, where many muon tracks 
occur. 

The above discussion deals with the ability to measure charge in an ideal 
detector. In a real detector, there are a number of reasons that we might 
assign an incorrect charge sign to a track. We separate these errors into 
two kinds which we call bias and randomization. One kind of error involves 
biases towards one sign or the other. Other errors will cause tracks to be 
randomly assigned as positive or negative. These randomized events will 
have an apparent charge ratio of unity. Since the true charge ratio is greater 
than one, any such errors will lower the measured charge ratio. We next 
describe two instances of each kind of error. 

The first kind of bias is associated with the magnetic geometry of the 
detector, which we call acceptance effects. Ignoring the very small magnetic 
field in the air outside the detector, the charge ratio of tracks which enter 
the detector is not affected by the magnetic field and/or the geometry of 
the detector. Instead, there are effects which have to do with cuts that are 
made to require a good track, and the efficiency of those cuts. From a given 
direction, tracks of one kind (the other) will be bent such that the track 
length in the detector is longer (shorter). In order to require a minimum 
number of hits to make a good charge determination, a minimum length cut 
may be made. Thus tracks in some directions may be biased towards one 
sign or the other. In the MINOS geometry, such directions occur as peaks 
or dips in the azimuthal distributions of the charge ratio. They tend to be 
symmetric with respect to the detector, but since a varying overburden can 
cause more tracks from some directions than others, they do not necessarily 
cancel. 

A second kind of bias comes from any misalignment of detector planes 
which causes a curvature in the coordinate system. A result of any such 
misalignment would be that straight tracks appear to have curvature. With 
the MINOS geometry, it was found that a systematic misalignment of even a 
fraction of a mm (as measured by the track's sagitta) would have a noticeable 
effect on the charge ratio. A curvature cut does not remove such a bias. 

Bias effects as described above are correctable in principle, with an ac- 
curate Monte Carlo. Acceptance issues can be corrected if the model of the 
magnetic field is accurate, and if the overburden is correctly modeled so that 
the expected fluxes as a function of zenith and azimuthal angle are known. 
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Misalignments can be modeled in principle, but must be done with a correct 
three-dimensional model of the detector, with an accuracy that is extremely 
difficult to achieve in practice. 

An easier way to deal with the above bias effects is to reverse the magnetic 
field. If the reversal of the field is done with sufficient accuracy, then all 
forms of bias described above cancel when one uses the geometric mean of 
the charge ratio for forward and reversed field running 10| 

[^forward ^ '^reversed] ^ (5) 

This cancellation is apparent when one considers the ratio of the positive 
muons in forward(reverse) running to negative muons in reverse(forward) 
running. The equality of these two values is a powerful consistency check. 

By contrast, effects which tend to randomize the charge ratio do not 
typically cancel in this manner and must be dealt with separately. An obvious 
type of randomization occurs for high momentum straight tracks above the 
MDM. If a track is straight within the resolution of the detector, then we 
will be unable to accurately determine the charge sign. For a straight enough 
track, a program will measure curvature either way an equal fraction of the 
time. The ffi that calculates the curvature can also calculate the error on 
the curvature. This form of randomization can be minimized by cutting 
on the fractional curvature error a{l/p)/{l/p) as is shown in Figure 5 of 
Reference However the fraction of mismeasured tracks may not follow 
a Gaussian distribution for a{l/p)/{l/p) and Monte Carlo estimates of the 
misidentification must be analyzed with care. 

A second form of randomization, which was observed in MINOS data, was 
generally associated with extra hits. These might be due to noise, crosstalk 
or demultiplexing issues, or also delta rays or other particles associated with 
the muon. For beam experiments, the magnetic field is perpendicular to the 
beam direction, to maximize bending power from the magnet. For cosmic 
rays, muons are observed at all angles. As described above, muon tracks will 
be bent much more in a toroidal magnet if they impinge on the detector with 
a low impact parameter. For those cosmic ray muons which are traveling 
mostly parallel to the magnetic field (at large impact parameter), an inter- 
esting situation can arise if the apparent bending due to multiple scattering 
is comparable to the bending from the magnetic field. If that happens, a 
reconstruction program which takes into account multiple scattering can get 
a good fit for a straight track for any momentum. The inclusion of a hit or 
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hits which do not belong on the track can then lead to a class of muon tracks 
which are really at high momentum, but get reconstructed at low momentum 
with a random charge sign depending on which side of the track the extra 
hits appeared. 

We can quantify this result as follows: Projected onto a plane, the RMS 
multiple Coulomb scattering angle can be written 

< Omcs >= [0.0136 /p{GeV/c)] * ^{L{meter)/X) (6) 

where X is the radiation length. The magnetic bending for this same length 
of track is given by 

OBend = L{meter) * BpERp{Tesla)/[3.33Q x p{GeV/c)] (7) 

Note that in the ratio OMcs/(^Bend, P cancels. For typical muon tracks in 
MINOS, the ratio rises from 0.14 at low impact parameter and ip to near 
unity at large impact parameter. Thus the type of randomization described 
in the previous paragraph can be preferentially removed by discarding events 
with large values of < 6mcs > /dsend- 

The cancellation of bias effects can be checked by studying distributions 
which have bias effects in them before and after the forward and reverse field 
data are combined. We have no such consistency check for randomization 
effects. Thus, the ability to detect and control all bias effects using the 
combined forward and reverse field data is stronger than the ability to cancel 
all randomization effects. This might lead to a tendency for experiments to 
report a charge ratio which is systematically smaller than the true value. We 
also note that not every experiment listed in Section [2] has taken data with 
a reversed magnetic field. 

4. Issues with Interpreting Data Obtained Underground 

4-1. The Relevance of Energy Loss Differences 

The charge ratio measured in underground experiments is a ratio of /i"*" 
to /i~ at the detector. For comparison with production models, the relevant 
quantity is the ratio in the upper atmosphere. If the energy loss of positive 
and negative muons were identical, as they go through the earth, these two 
ratios would be the same. The leading energy loss processes of and /i~ 
are the same of course, but there are corrections of the order of the fine 
structure constant a. We evaluate these small differences, first for ionization 
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energy loss, and then for radiative processes. The latter difference is found 
to be negligible. The small difference in ionization energy loss is used to 
calculate the difference in range and the correction to the charge asymmetry 
as a function of slant depth in Section 14.41 

The statistical energy loss of muons, traversing an amount X of mat- 
ter in g/cm^, with energies far above the Bethe-Bloch minimum is usually 
parameterized as 



dX 



(8) 



where a is the collisional term (i.e. ionization) and h in the second term ac- 
counts for the three radiative muon energy loss processes: Bremsstrahlung, 

these 



pair production and photo-nuclear interactions. In Table [Tjl8|[19 
energy loss parameters are listed for standard rock. The critical energy 
where ionization losses equal radiative losses in standard rock is approxi- 
mately 0.6 TeV. The average muon surface energy for a muon which reaches 
2000 mwe is greater than 1 TeV, so the h term and its energy dependence are 
important in calculating the energy loss. We have investigated differences in 
the a and h terms for /i"*" and ^~ . 
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Table 1: Average muon energy loss parameters calculated for standard rock [18|] 



4-2. Difference in Ionization dE/dx for and fi 

At low energies, around the Bethe-Bloch minimum, the difference in ion- 
ization energy loss is known as the Barkas effect 2C [, and there have been 



efforts to both measure and calculate the difference [2l[|. Calculations show 
that negative particles lose energy at a slower rate, with the difference drop- 
ping from tens of percent at MeV energies to about 0.3 % in the GeV range. 
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Figure 2: Calculated difference in ionization energy loss between positive and negative 
niuons in standard rock (average nuclear properties: Z — 11, A = 22 [Toj). 



Such differences were experimentally verified both at MeV energies 22| 23|] 24 



and in the GeV range [25[. At higher energies, this difference in ionization 
energy loss has usually been neglected, and we are not aware of any measure- 
ments. As described in Reference j26[, the usual ionization energy loss term 
for muons (of either sign) depends on 2;^, and the difference between /z"*" and 
/i~ arises from a small additional correction term. This correction term 
in dE/dx is: 

\UeV,.a',r^^ (9) 

where a is the fine structure constant, z is the charge, /? is the relativistic 
velocity, and Z and A are the nuclear properties of the material through 
which the muon is passing. The absolute value of the difference in ionization 
energy loss between positive and negative muons in standard rockflSJ] is plot- 
ted in Figure [2J It is fairly constant above 10 GeV, at a value corresponding 
to approximately 0.15 % of the mean energy loss in the ionization dominated 
energy regime (c./. Table [T|) (irf . 

Calculated Difference in Bremsstrahlung dE/dX for /i+ and fi~ 
Above an energy near 0.6 TeV in standard rock, radiative energy loss be- 
comes comparable to ionization energy loss, and continues to grow at higher 
muon energies. Radiative energy loss has been calculated using the Eikonal 



approximation in References [28| and [29| . For the difference between positive 
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Figure 3: Calculated fractional difference in Bremsstrahlung energy loss between positive 
and negative muons. 



and negative particles, the leading term cancels, and only a term proportional 
to M/E = 1/7 survives, which does depend on the sign of the charge 3ol |: 



rdE_-ifi+ _ ldE_lfJ'- R7rv 

\-dX^br(i'ms IdX^brems OZjUl ("SfW 

where 7 is the Lorentz factor of the muon. Again, the /x"*" has a slightly higher 
energy loss. This fractional difference decreases with energy and is already 
negligible where radiative energy losses become important. This fractional 
difference is plotted in Figure [3] for muons in standard rock. Presumably, 
the same fractional difference can also be assigned for pair-production, as 
the underlying process is a two-photon exchange between the muon and the 
constituents of the nucleus, and thus the cross sections for /x"*" and /i~ should 
scale in the same way as for Bremsstrahlung. 

4.4- Range Underground and Muon Charge Asymmetry 

Taking the vertical muon intensity from the Gaisser parameterization 
of the muon flux at the surface (see Equation [2]) and propagating this en- 
ergy spectrum underground according to statistical ionization and radiative 
energy losses, one can calculate the underground muon intensity. This pro- 



cedure is described in detail in [31I] for overburdens of standard and Soudan 
rock (MINOS). First, the average muon range underground, for each value of 
surface energy, is computed. For this, the energy dependent a and Th values 
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were parameterized for standard and Soudan rock as in [31[. The additional 
ionization loss according to Equation M was added for /i"*" to the value of 
the function for a (subtracted for Radiative losses (90% of S6) were 

scaled with the energy dependent fractional difference SZa/j according to 
Equation [To] as = 0.9 ■ S6- (1 ±4Za/7) for /x"*" and respectively. The 
contribution from photo-nuclear production (DIS), which we are not aware is 
charge dependent, was not scaled. For each initial value of muon energy, the 
slant depth in meter-water-equivalent where the muons of different charge 
range out was determined. 

We computed the underground intensities of positive and negative muons 
as a function of slant depth for a given rock composition. This ratio of the 
/i"*" and /i~ intensity curves is shown in Figure H] for Soudan rock. The upper 
curve corresponds to the fractional difference in integral intensities of /i"*" 
and fi~ at a given slant depth. For slant depth values above about 1000 
mwe the underground ratio N{^^)/N{fi~) is lowered by roughly 0.4%. The 
detected intensity corresponds to the charge ratio of the muons at depth 
below the MDM for MINOS. For increasing slant depth values the measured 
underground ratio N{fi~^)/N{fi~) is further reduced and saturates at about 
0.6 % below its surface value for slant depths larger than roughly 5000 mwe. 

This value can be qualitatively understood as follows. If the muon flux 
was proportional to and if the muon range was strictly proportional 

to energy, then less positive muons would survive underground and 

^underground ^surface ^'^l^ ^(-^m)) fill 

So the change in the charge ratio is A(r^) ~ 3.7AEfj^/E^ Thus, the small 
fractional difference in energy loss for /z"*" and /i~ of the order of 0.15% 
predicted by theoretical calculations at high energies, gets amplified by a 
factor of about 3.7. The impact of the rock composition is almost negligible, 
as the induced muon charge asymmetry under Soudan rock lowers the surface 
value of the ratio by an additional amount less than 0.02 % compared to 
standard rock. Figure H] can be used to correct the underground measured 
muon charge ratio to its surface value. For MINOS, this corresponds to a 
correction which increases the measured value by 0.6%. 
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Figure 4: Calculated ratio of positive to negative vertical muon intensities in Soudan rock 
as a function of slant depth for equal surface intensities. The upper curve is for all muons, 
the lower curve is for muons with a remnant momentum of less than 250 GeV/c the 
maximum detectable momentum in the MINOS Far Detector). 

5. Model of Pion and Kaon Contributions to the Charge Ratio 

We have investigated a generalization of Gaisser's Equation |2] to study 
separately the positive and negative muon intensities. Using the positive 
fraction parameters f-„ and fx, the energy dependency of the positive and 
negative muons is given by 



where and ex have been replaced by their numerical values. We can use 
equations [12] and [13] to calculate the surface muon charge ratio: 



The charge ratio of muons from pion decay is r^^ = fn/{i — fn) and 
from kaon decay is = — fx)- We will refer to the implications of 





(13) 
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Equation [T¥l with energy independent parameters as the "pika" modeL There 
are several interesting features of this modeh 

1. The relative intensity of cosmic ray pions and kaons that contribute 
to muon production can be extracted from surface and underground 
muon charge ratio experiments. 

2. The muon charge ratio does not depend upon the muon energy and the 
zenith angle separately, but on the product E^^^"-'^'^ cos 9. This product 
of terms controls the relative portions of interaction and decay for both 
pions and kaons. At a fixed value of E^^f-^ cos 9 , the intensity ratio 
of muons from pions and kaons is constant. 

3. There are contributions to the muon charge ratio from both pions 
and kaons from 1 GeV to 10 TeV. The contribution to the ratio from 
kaons does not vanish at GeV energies - it is just smaller than at TeV 
energies. This effect can be seen in Figure |5] which displays the contri- 
bution to the muon flux from each of the four mesons. 

4. The pika model postulates an energy independent r^r = tt^/tt" ratio 
related to and energy independent = / K~ ratio related to 
/x, and an energy independent ■n/K ratio embodied in the Gaisser 
constant 0.054. 

5. Contributions from charm particle production have been ignored be- 
cause their effect on the ratio below 10 TeV is expected to be negligible. 

Where Feynman scaling is valid, the fraction x = Emeson/ Eproton for vr 
and K secondaries does not depend upon Eproton- Then and fx are also 
energy independent. Here we explore features of our simple model, and in 
the next section we shall compare the qualitative features of our model with 
full simulations of the charge ratio. 

Our parameterization yields additional insight into the sensitivity of var- 
ious underground experiments which measure the muon charge ratio. Under- 
ground experiments can effectively observe a much smaller range of E^^^"''^'^ cos 9 
than a range of E. To see this, consider a simple approximation with five as- 
sumptions: 

1. The parameter a in Equation [8] is constant 

2. The parameter b in Equation [8] is constant and set to zero. 

3. The earth/rock density above the underground detector is constant 

4. The surface above the detector is fiat 



15 



0.6 

c 

§ 0.5 
E 

1 0-4 
o 

E 

2 0.3 



o 

■■s 0.2 



i 0.1 

o 




3C+ muon flux contribution 
n- muon flux contribution 
K+ muon flux contribution 
K- muon flux contribution 



-I I I I I I I I I I 



10"^ 10"' 



10"^ 1 
E^,cos(e) [TeV] 



10 10' 



Figure 5; Contributions to the muon flux from the four mesons using the pika model 
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5. The MDM in the underground detector is neghgible compared to the 
energy loss from the surface to the detector. 

The first four assumptions are equivalent to stating that the muon energy 
loss in the earth is proportional to the overburden which only depends upon 
the muon zenith angle. In that case the energy loss through the overburden 
is Eioss = Emin/cos{d), where Emin is the minimum energy loss for a vertical 
cosmic muon. The fifth assumption implies that the underground energy of 
the muons used for the charge ratio measurement is much smaller than their 
surface momentum. 

With these five assumptions, the muon intensity distribution in E^^^^"''^'^ cos 6 
would be measured to be a delta function at a value of -Emm- While at large 
zenith angles the surface energy increases due to the l/cos{9) dependence of 
the slant depth, the combined quantity E^^'^^'^'^^ cosO remains constant. It 
is illuminating to compare this naive prediction to the actual distribution 
in the MINOS detectors when no such assumptions are made. These are 
shown in Figures E] and [71 Both of these MINOS E^^^^"''^'^ cos6 distributions 
are considerably narrower than the corresponding Egurface distributions. The 
b(E) radiative term yields the largest contribution to the width of the mea- 
sured Efl^^f"''^^ cos6 distribution for the MINOS Far Detector. In the Near 
Detector distribution, the largest contribution is the larger ratio of maximum 
detectable momentum to energy loss in the overburden. 
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Figure 6: Distribution of E and 
^surface q f^^, Mj^QS data muons in 
the Far Detector, after cuts (GeV). 



Figure 7: Distribution of E and 
^surface g, f^j. ^INOS muons in the 
Near Detector, after cuts (GeV). 



We have used Equation [H] and the measured muon charge ratio to study 
and tk- We have done chi-squared fits in E^^^^"''^^ cos9 to the MINOS 
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Figure 8: Pika model fitted to L3+C and MINOS (Near and Far) data sets 



Near Detector and Far Detector data, and to the L3+C data. The fit yields 
U = 0.5510 ± 0.0006 and fx = 0.7006 ± 0.0061. These values lead to a 
muon charge ratio from pion decay of r^^ = 1.227 ± 0.003 and a muon charge 
ratio from charged kaon decay of tk = 2.34 ± 0.07. The errors in the two 
parameters are highly correlated. This fit is shown in Figure E] along with 
the MINOS (Near and Far) and L3+C data. The agreement between the 
parametrization and the data is excellent. 

Using the f values from this fit, we plot Equation [H] for an extended 
range of muon energy in Figure M The critical energies discussed above for 
pion and kaons are indicated by the arrows. The low energy asymptotic 
ratio is (/^ + ?7//^)/[l — + ri{l — f^)] = 1.26. The high energy ratio is 
(ctt/tt + ei^^?//^)/[e7r(l - /tt) + e/^^?(l - fx)] = 1-45. The low and high energy 
asymptotic values both include muons from vr decay and from K decay. 

The fit was repeated using just the smaller data sample of the MINOS 
Near DetectorQ and Far Detector. We obtain = 0.5538 ± 0.0070 and 
/k = 0.693 ± 0.027, which imply = 1.241 ± 0.035 and = 2.26 ± 0.29. 

The parameterization we have used seems sufficient to represent the pub- 
lished data sets. 

Our fits to r,r gives values near expectations. Our fit to tk = K^/K~ in 
atmospheric showers yield values near 2.3. It is clearly difficult to directly 
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Figure 9: The pika model's muon charge ratio for an extended range of energy. The critical 
energies for pion and kaon decay vs interaction are shown by the arrows. 



measure the atmospheric kaon charge ratio. Equation [H] provides a well 
defined parametrization for future studies of this subject. 

6. Comparison with Full Simulations and Calculations of Atmo- 
spheric Muons 

Full simulations of the muon charge ratio have been published by Honda 
34 1 and CORT j^; a full analytical calculation of the ratio has been published 



by Lipari 3^. These model results of the ratio versus the muon surface energy 
j^surface displayed in Figures [TOl [11], and [12] for various ranges of cos(6'). 
The Honda results show small dependencies upon the zenith angle, and an 
interesting dip in the ratio near 200 GeV. The CORT results show a larger 
dependency upon the zenith angle as the energy increases. The Lipari results 
show even a larger dependency upon zenith angle at energies above 10 GeV. 
The results from the models differ noticeably from each other and from the 
MINOS and L3+C data sets. 
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Figure 10: Honda calculations of the muon charge ratio at 5 zenith angles 
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Figure 11: CORT calculations of the muon charge ratio at 4 zenith angles 
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Figure 12; Lipari calculations of the inuon charge ratio at 5 zenith angles 
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In our analysis, we do investigate (1) the dependency of these charge 
ratio simulations/calculation on the E'l^'^^"''^'^ cos9 variable discussed in the 
previous section of this paper, and (2) their consistency with the pika model 
of the charge ratio. Our observation will be that E^^^^"''^'^ cos6 is a more 
useful variable for ratio analysis than just the muon energy E^^^-^"''^'^. 

Figure [12] displays Honda's muon charge ratio simulations as a function 
of E^'^-l'"-'^'^ cos6 using published values from five ranges of cos(^) from 0.3 to 
1.0. For a wide range of the variable E^^^^"''^'^ cos9, the charge ratio is nearly 
independent of cos(^), consistent with the representation of the pika model. 

Figure [H] displays the corresponding CORT theory simulations of the 
charge ratio as a function of i^^"''/'^^^ cos 6 using four published values of 
cos(^). In the E^^^^"''^'^ cos6 range from 5 GeV to 600 GeV, the simulations 
show minimal dependency on cos(6'). Above 600 GeV, the charge ratio model 
results do increase with increasing values of cos(^). The prediction of the pika 
model that the ratio is independent of E^^^^^"-'^'^ cos 6 is consistent with the 
CORT simulation for a remarkable range of E^^^^^"-'^'^ cos 6. 

Figure dS] displays the Lipari calculation of the charge ratio as a function 
of £;^"''/'*ce cos^ for five values of cos(6') from 0.2 to 1.0. It is apparent that a 
significant portion of the ratio variation with cos(6') is minimal when the ratio 
is displayed as a function of i^^"''/'^'=e (^Qg^_ Again, the pika model prediction 
of a charge ratio independent of E^^'^^'^'^^ cos 9 is consistent with the Lipari 
calculation over four orders of magnitude in i<;^«^/'ice q This agreement is 
to be expected and is due to Lipari explicitly using the identical E^^'^"-'^^ cos 9 
dependency in his muon flux calculation, as Gaisser does in his model ^ (our 
Equation [2]). 

The observed dependence of the above three models on £;«"»"/<^ce ^^^g q 
vides support for analyzing the charge ratio data using the pika formula. 
Of course it would be interesting to understand the source of the remaining 
variations in the simulations/calculation as a function of E^^'^-^"''^'^ cos9. 

Next we quantitatively compare the pika formula to the three simula- 
tions/calculation discussed above. Two parameter chi-squared fits were per- 
formed to the simulation's/calculation ratio values assuming equal uncer- 
tainties on each point within a model. The best fit for each model is shown 
as a smooth line in the corresponding figure. For the CORT model, the 
pika model fits well the simulation from 20 GeV to 600 GeV at all angles. 
The fit yields f-„ = 0.557 and fx = 0.705. For the Honda simulation, the 
pika model can not reproduce the simulation's dip at ~200 GeV, although 
it does describe the simulation's overall energy dependency from 5 GeV to 
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5 TeV. The Honda fit yields /^r = 0.5615 and fx = 0.6207. For the Lipari 
calculation, the pika model must of course reproduce well the full calculation 
from 3 GeV to 10 TeV due to Lipari's usage of E^^'^^'^'^^ cos 6. The fit yields 
/tt = 0.5551 and fx = 0.7413. To summarize, we conclude there is physics 
content in using the variable E^^^^^"-'^^ cos 6 in studying the results of these 
three simulations/calculations. However, the pika formula can not reproduce 
the full richness of these simulations, and of course the fitted pika parameters 
are significantly different among the three simulations. 

To complete this study, we have compared the L3+C and MINOS charge 
ratio data sets to the same three simulations/calculation. So as to concisely 
represent the three simulations/calculation, we have used the fits of the three 
models to the pika formula to compare with data. As seen in Figure dni none 
of the models (as parameterized by the pika formula with varying degrees of 
success) agree with the high E^^^^"''^'^ cos 6 cosmic ray data. Comparing the 
pika model results shown in Figure [9] and Figure [161 it is apparent that none 
of the three simulations/calculation are using optimal values of /^r and fx- 
This can be seen in Table E] and is discussed further in Section [Si 



7. Effect of Helium on the Charge Ratio 

Another process that could affect the energy dependence of would 
be a different spectral index for the heavier cosmic ray primary flux than for 
the proton flux. This would introduce an energy dependence to the incoming 
proton to neutron ratio. It is not well established that the energy dependence 
of the heavy primary intensities is different than that for Hydrogen in the 



10 TeV energy range. For this study, we will use Gaisser and Honda's [35 
parameterizations of the primary flux as a function primary component k 
with energy E^, given as 

(f){Ek) = K{Ek/GeV + b x exp[-c^/ (Ek/GeV)])-" (15) 

The parameters in the above equation are given in the Table [2] below. Note 
that this reference has a spectral index decreasing slightly with increasing 
primary mass. 

Figure [T7| displays the contributions to the flux as a function of the kinetic 
energy per nucleon. Based on Monte Carlo calculations using CORSIKA, 
the mean surface energy of muons which reach 2100 mwe underground is 9% 
of the primary nucleon energy. The mean fraction of the primary Helium 



24 



1.6 
1.55 

1.5 
1.45 

1.4 
1.35 

1.3 
1.25 

1.2 



□ cos e = 0.5-0.6 
A cos e = 0.6-0.7 
i COS e = 0.7-0.8 
« COS e = 0.8-0.9 
* cose = 0.9-1.0 
Honda fit: fpi=0.561 5 and fka=0.6207 




_i I I I I I I I I I I I 



10' 



10"^ 10"' 

E^cos(e) [TeV] 



10 



Figure 13: Honda calculations of the muon charge ratio vs £:s"'"/«ce ^ 5 zenith angles. 
The curve is a fit of these simulation results to the pika model. 
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Figure 14: CORT calculations of the muon charge ratio vs _E«"'"/«'=e ^ ^t 4 zenith angles. 
The curve is a fit of these simulation results to the pika model. 




26 



1.6 
1.55 

1.5 
1.45 

1.4 
1.35 

1.3 
1.25 



cos=1 .00 
cos=.60 
cos=.40 
cos=.30 

cos=.20 

Lipari fpi=0.5551 and fka=0.741 3 



•'A 



o. 



-o ■ A 



1.2^ 
10' 



_i I I 



_i I I 



_i I I 



-I I I 



10"^ 10"^ 1 

E^cos(e) [TeV] 



10 



Figure 15: Lipari calculations of the muon charge ratio vs E^'^^""' cos ^ at 5 zenith angles. 
The curve is a fit of this calculation's results to the pika model. 
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Figure 16: Fits of the pika formula to the three models of the charge ratio vs cos 9 

compared to L3+C and MINOS data sets. 



Parameter 


a 


K{m?' s sr GeV) ^ 


h 


c 


Hydrogen 


2.74±0.01 


14900±600 


2.15 


0.21 


He (A=4) 


2.64±0.01 


600±30 


1.25 


0.14 


CNO (A=14) 


2.60±0.07 


33.2±5 


0.97 


0.01 



Table 2: Primary flux parameters used in the text. 

energy, per nucleon, that is transferee! to the muon is 11%. We will assume 
the fraction of the energy transfer for primary carbon, nitrogen, and oxygen 
to be identical to that of Helium. In the following calculations, the energy 
profiles of the protons and Helium from Monte Carlo are used instead of 
these average values, although the effect is not large. 

We will calculate the corrections to the fitted and fx parameters in 
Equation [T3] to account for the possible energy dependency of the incoming 
proton to neutron ratio (due to the Helium and CNO flux). We use a model 
to estimate the impact of both pion and kaon production by protons and 
neutrons. We assume two contributions for meson production: pair produc- 
tion of equal number of positive and negative mesons, and leading particle 
production using quark counting to estimate the excess. This calculation will 



28 



Hydrogen 

Helium 

CNO 



-I I 



_i I 



_i I 



_i I 



_i I 



10' 



10"^ 10"' 1 10 

kinetic energy per nucleon [TeV] 



10^ 



Figure 17: Fractional contribution of nuclei to the all nucleon spectrum 
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use meson charge ratios near the previous fit values. We start with = 1.25 
and tk = 2.61 and see how much they change due to the different spectral 
index for H and He. 

First we will model the change in the muon charge ratio due to just pion 
decay when Helium and CNO are added to the primary proton fiux. Using 
CORSIKA separately for cosmic ray protons and Helium of energy Ecr, we 
have calculated weights wMn 11 bins of the ratio E^/ Ecr. In a simple quark 
model the leading proton and neutron particle in the forward direction will 
give rise to a charge ratio of 1.00 for symmetric nuclei A-Z=Z. 

Setting the muon charge ratio equal to 1.25 at 1 GeV and fixing the 
normalization of the Hydrogen fiux to be 1.317 (in order to get back the 
charge ratio of 1.25 at 1 GeV), one has for the muon charge ratio 

10 

r-(E^) = 5^[0.5x (ti;; + <Jx 

i=0 

Wil.m X 0p(Ecfi) + <,1.0 X A<j)He{EcR) + <el-0 >< M^CTVO (^Cfl) 

+ W^HMueiEcR) + W^Hel^CNoiEcR) 

(16) 

The 1.0 values are due to the equal production rates of leading tt"*" and 7r~ 
in our quark counting model. 

A subtlety of our notation is that the vr (and K) superscript represents 
the energy dependent correction to the charge ratio arising from the different 
possible spectral index of heavy cosmic ray primaries, while the subscript 
represents the assumed energy independent meson ratio. Figure [TS] displays 
the corresponding fractional change in the charge ratio as a function of the 
muon surface energy. Note that the effect is about a 3 percent reduction at 
1 TeV. This dependency can be fit to a polynomial in log(E), also shown in 
Figure [T8| with the result 

r^(E^)/1.25 = 1 - 0.00575x/o^io(E^/GeV) - 0.00155x/o/io(E^/Ger) 

(17) 

Next we include the contribution from kaon production by the primaries. 
The cross section for associated production of KK^ is much larger than that 
for K^K~ pairs for both incident protons and neutrons. Again assuming 
a quark model for the leading kaon in the forward direction from incident 
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protons and neutrons, one only gets leading , but twice as many from 
protons than from neutrons corresponding to their number of u quarks. 

p{uud) + air^K^ {us) + ■ ■ ■ (18) 

and 

n{ddu) + air-^K^{us) + ■ ■ ■ (19) 

In this model, the positive kaon excess originating from protons {r^ — 1) is 
twice as large as the excess originating from neutrons. For symmetric nuclei 
with A-Z=Z this leads to a charge ratio 

1 + i^^(r^ - 1) = 1/4 + 3r;,/4 (20) 

We further assume that the muon charge ratio at 1 TeV is 1.374 (near the 
MINOS value). This then yields 

10 

r^(E^) = ^[0.5x(^; + <Jx 
wp.l X (t)p{EcR) + <e2-275 X A(j)He{EcR) + <e2-275 x IAcPcno{Ecr) 

wi4)p{EcR) + w\j^A(t)He{EcR) + w\j^U(j)cNo{EcR) 

(21) 

The value 2.275 is due to Equation [201 The leading particle positive excess 
dominates over K^K~ pair production, unlike in the pion case. Again we 
parameterize the dependency in IoqiqE^ from Equation [17] such that 

(1 - r^(E^)/2.61) = 2/3 X (1 - r^(E^)/1.25) (22) 

The effect on the muon charge ratio due to the heavy ions requires that 
the pion and kaon fractions be modified as follows: 

f:iE,) = l/il + l/r^xr-iE,)) (23) 

r^{E,) = 1/(1 + l/rKXr^{E,)) (24) 

Including these contributions from heavy primaries does have an impact 
on the calculation of parameters of the pika model. The charge ratio no 
longer depends just upon E^cosO. While this effect is small it is certainly 
present and can be accounted for in modeling. Figure [TH] shows the effect 
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Figure 19: Fit to the meson charge ratios modifying the pika model with the energy 
dependence of Helium and CNO production. 



of the heavy primaries on the two parameters of the pika modeL The fit 
to the data is almost indistinguishable from the previous fit. Note that an 
increasing heavy primary fraction at high energy will decrease the charge 
ratio, so to fit to the high MINOS points, a larger value of tk is required. It 
is clear that simulations of the muon charge ratio need to include the possible 
different energy dependence of heavy primaries. 

We repeat that the choice of spectral index in Table [2] is to illustrate the 
size of a possible effect. If it turns out that the spectral index is independent 
of chemical composition, the effect described in this section will not exist. 

8. Discussion 

We have considered the change in the charge ratio from the fact that 
fi~^ lose slightly more energy than fi~ while penetrating the overburden of 
an underground detector, both for ionization and catastrophic energy loss. 
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For catastrophic energy loss, the effect is neghgible, even for the small errors 
in MINOS. For ionization, we have calculated that the 0.15% difference in 
energy loss leads to a 0.6% difference in measured at 2000 mwe. This 



effect was mentioned in Reference [1| but not used to correct the reported 
charge ratio there. In this paper, we correct the MINOS data to slightly 
higher values. For the Far Detector data reported in Reference this is a 
+ 0.6% correction, and for the near detector at 400 mwe the correction is 
+0.37%. We estimate a +0.29% correction for L3+C. 

The central result of this paper is the development and application of 
Equat ion [T^ and the consequent relationship of to and tk- We have used 
that equation to study its consistency with several simulations/calculations, 
and to make fits to data. In Table [21 we show the fits for our parameters to 
the three more detailed models. We also show the fits to the data reported in 

and 36 1 . One can also compare these values to those 
and Agrawal ^\ ; assuming a proton primary cosmic ray 



References |2[[ 
used by Gaisser 

flux (ignoring the Helium component), these models use published pion and 
kaon production data from accelerator proton beams to specify the and 
parameters. In particular, Gaisser specifies = 1.4 and =3.2; Agrawal 
specifies = 1.35 and tk =2.92. 







tk 




Calculation 


Lipari 


2.87 


1.25 


Simulation 


CORT 


2.39 


1.26 


Simulation 


Honda 


1.63 


1.28 


Data fit 


MINOS and L3 


2.34 ± 0.07 


1.227 ± 0.003 


Data fit 


MINOS N+F 


2.26 ± 0.29 


1.241 ± 0.035 


Data fit 


Hehum 7 2.74 ^ 2.64 


2.73 ± 0.09 


1.234 ± 0.003 



Table 3: Pika fits to tk and for a calculation, simulations, and data. 

We have shown the importance of using £;^"''/'^ce ^-.^g q (^Y^stea.d of just 
energy) in the analysis of the atmospheric muon charge ratio for energies 
above 10 GeV. Understanding of the muon charge ratio would benefit if future 
experiments and theory simulations would provide results using £'^"^/'^ce q 
The rise in the charge ratio, which occurs between the vr and K critical 
energies of 115 GeV and 850 GeV, can be used to determine the meson 
charge ratios /K" and jn^ . This is an important new method for 
obtaining information on these ratios. 



34 



We further showed the effect on the charge ratio assuming that there 
is a spectral index for cosmic ray Hehum nuclei which is different than the 
spectral index for protons. Using an index of 2.74 (2.64) for H (He), we get 
a change in tk from 2.3 to 2.7 compared to an analysis where the spectral 
indices were the same. 

Both the data and the simulations seem to have consistent values for r^. 
A general feature of our analysis is that the new MINOS data by themselves 
suggest a lower value of rx than is used in many of the simulations. How- 
ever, it is important to note that the systematic effects considered in this 
paper, which include remaining errors from randomization, corrections for 
differences in dE/dx for /i+ and /i~ and a possible lower spectral index for 
Helium, would all tend to raise the fitted values of rx- 

Several effects which have not been explicitly considered here are ex- 
pected to be small. These include the production of muons from charm and 
other heavy particles, components of the cosmic rays heavier than Helium 
and possible differences in their spectra, and a variety of scaling violations 
which would have the effect of making f^^ and fx energy dependent. These 
effects need full simulations to evaluate fully; but it is important that the 
full simulations yield the experimentally measured average values of these 
parameters. In that context, the analysis presented in this paper can be 
useful. 

9. Summary 

We have reviewed several factors which affect the measurement and in- 
terpretation of the muon charge ratio deep underground. While many 
quantities in cosmic ray physics are difficult to measure precisely, has 
been measured in MINOS with a statistical accuracy better than 0.3%. In 
order to take advantage of this high precision, an experiment must control 
systematic errors to a comparable level, which is challenging. We pointed 
out two kinds of systematic error, those which might bias the measurement 
of r^, and those which randomize (and hence lower the measured value from 
the true value). Bias errors can be canceled to high precision by using data 
with both magnetic field polarities and Equation and the success of this 
cancellation can be checked for consistency. Randomization errors, on the 
other hand, require severe cuts that affect the statistical precision. Variables 
which may be related to randomization effects can be identified and used 
to cut out data which may be affected by randomization. However, there 
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is no independent way to determine if all such effects have been eliminated. 
The possibility exists that the true values of the charge ratio are higher than 
those that have been reported. 

The MDM is an important parameter for an experiment with a magnet. 
It affects both the calculation of the systematic error and the energy resolu- 
tion, and also directly limits the range of detected muon energies useful for 
measuring and £'««»'/ace 

We note that to fully explore the rise in the charge ratio observed by 
MINOS, there is a need for additional precise ratio data in the E^J^'^^"''^'^ cos9 
ranges of 0.2-0.8 TeV, and above 3 TeV. 
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